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We use the simple instanton motivated Nambu Jona-Lasinio - type model to calculate a twist 
3 pseudo-scalar pion light cone wave function ip^ s . Using the normalization condition for ip^ s we 
calculate the quark condensate and also the gluon condensate, which agree with the phenomenolog- 
ical values for these quantities. Since we can compute also the k% dependence of the light cone wave 
functions, we calculate k\ moments of the pseudo-scalar and axial-vector wave functions which are 
related to the mixed vacuum condensates. This allows us to extract the condensates and compare 
, them with existing estimates. 
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I. INTRODUCTION 

, Soft pion theorems provide link between dynamical objects like the light cone wave functions fl]-^| and static 
' properties of the physical vacuum MM (for review see §). Customarily in QCD pion properties are expressed in 
^sO ■ terms of the nonperturbative condensates, which are estimated by some theoretical model, like the QCD sum rules 
0^ ' H (see for example Q), or the instanton model of the QCD vacuum pO| , pl| (for recent review see e.g. [^2|). Here we 
propose to invert this logic. In Ref. we have discussed a simple and tractable model, which describes very well the 
leading twist pion wave function V'^not only in longitudinal variable it, but also in kx- In this paper we calculate the 
twist 3 pion wave function ifj^ s pjl4|-[lq| and try to extract from both tjj A v and ip^ s the quark condensate and mixed 
condensates of dimension 5 and 7, and the gluon condensate as well. The model, which we use is a simple nonlocal 
chiral pion-quark model ]l7| , p^] , which in principle can be obtained from the instanton model of the QCD vacuum 
[ 12 1 . However, since we want to work directly in the Minkowski space, we have to modify the form of the nonlocality 
[13, p7|p8| in order to make the calculations feasible. The price is that, of course, we loose the direct contact with the 
Oh, original model. Therefore our aim is twofold. Apart from giving estimates of the condensates we want to perform 
tests in order to gain confidence in the model as well as to find its limitations. 

In most applications of the effective quark-meson theory described above, an approximation of the constant con- 
stituent quark mass was used. This approach has an impressive success in describing baryon (which emerge as solitons) 
properties (see |jl2| and references therein). Recently an attempt to include the momentum dependence in calculating 
solitons in the nonlocal NJL model has been reported in Refs. ]l9[ . 

The instanton model of the QCD vacuum is formulated in the Euclidean space-time, whereas the light cone wave 
functions are naturally defined in the Minkowski space-time. One can in principle Wick rotate the integrals and 
perform the calculations in the Euclidean metric, as has been for example done in Refs. (20 22 (for early attempts 
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to calculate the pion wave function in the instanton model see Ref. ]23|]), however in this work we shall perform all 
calculations directly in the Minkowski space-time. 

The instantons induce both gluon condensate, which is proportional to the density of instantons, and quark con- 
densate, which emerges due to the derealization of the quark zero modes Therefore both quark and gluon 
condensation occurs at the same scale Qo, which is associated with the average instanton size l/~p = 600 MeV. After 
integrating out gluons and performing the bosonization, the model reduces to the simple Nambu- Jona Lasinio type 
model where the quarks interact nonlocally with an external meson field U fl2| , p4| : 

si = J -j^mvw^u^(k-i)y/wm) « 

and C/ 75 can be expanded in terms of the pion fields: 

^ 75 = l + f l 5 r A n A -^ A n A + ... (2) 
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Here F n = 93 MeV and M(k) is a momentum dependent constituent quark mass. In principle M(k) has been calculated 
in the instanton model in Euclidean space time. Here, following Refs. ^^H , we wish to perform calculations directly 
in the Minkowski space. To this end we shall choose a simple pole formula |b|] 

/ A 2 X 2n 

which for n ~ 2 — 3 and for k 2 < reproduces the k dependence obtained from the instantons reasonably well Ei] . In 
order to check model sensitivity to the specific form of the cutoff function (||) we shall vary M from 325 to 400 MeV 
and n from 1 to 5. Notice that M(k) of Eq.(||) provides an UV cutoff for the loop integrals entering the expressions 
for the pion wave functions. This is true both in Euclidean and in Minkowski space. 

Using this simple prescription we study various condensates, which are related to the normalization and k\ mo- 
ments of leading twist axial-vector (AV) and nonleading twist pseudo-scalar (PS) pion light cone wave functions (for 



definitions see e.g. Ref. [[15 16 1) as well as the PS wave function itself. 

^From the point of view of QCD the quantities we calculate depend on a nonperturbative scale Qo which, however, 
must not be confused neither with the constituent mass M nor with an auxiliary parameter A, which is fixed by the 
normalization condition for ip£ v . For k 2 < the Ansatz (||) should imitate M(k) obtained from the instantons. And 
for the latter, as explained above, Qo ~ 1/p = 600 MeV. It is therefore natural to assume that Qo is of the order of 
a few hundred MeV irrespectively of M and A. The precise definition of Qo is only possible within QCD and in all 
effective models one can use only qualitative order of magnitude arguments to estimate Qo- A more practical way to 



determine Qo will be discussed in Sect.VC where we will associate Qo with transverse integration cutoff Kt which, 
as we shall see, is of the order of 760 < Kt < 1100 MeV. 

As a first step we calculate the quark condensate, which enters the normalization of the PS pion light cone wave 
function. The numerical result depends upon M and n and varies between —(258 Mev) 3 and —(328 MeV) 3 . The 
commonly used phenomenological value of (qq) is lower and reads approximately —(250 MeV) 3 . Therefore our highest 
value for the quark condensate overshoots the phenomenological value by a factor of 2. The reason for this rather too 
high values of (qq) is a relatively slow convergence of the dk T integral for this quantity. This is even more visible in 
the case of the k T moments of the PS wave function. Indeed, one can relate both (k T ) Av and (k T ) ps to the same 
mixed condensate of dimension 5, (igqa ■ Gq), which yields (k T ) Av / (k T ) p = 5/9 ||. In our case this ratio is two 
times smaller. It is, however, interesting to observe, that if we cut the dk T at some (fc^) max = K T , which we choose 
in such a way that (qq) — —(250 MeV) 3 , then the ratio of (k T ) is stable and approximately equal to 5/9. The mixed 
condensate (ig qa ■ G q) comes out between —(427 MeV) 5 and —(443 MeV) 5 in a fair agreement with the estimates 
from the QCD sum rules Ref. p5| . Unfortunately the higher moments of the PS wave function are still too large even 
with a k T cutoff. 

One of the advantages of our method is that the analytical expression for the quark condensate is given in terms 
of a Minkowskian integral, which in a limit of a constant M(k) and k 2 M — ► —k 2 E reduces to the well known Euclidean 
form. By comparing the two expressions one can by inspection guess a continuation prescription, which allows to 
rewrite certain Euclidean integrals as the Minkowskian ones. We use this, in some respect ad hoc, prescription to 
calculate the gluon condensate with very encouraging result: (a/wGG) = (378 MeV) 4 to (383 MeV) 4 . 

Finally we plot the PS pion light cone wave function for our set of model parameters. We find that <f)^ s (u) vanishes 
at the end points u = and u = 1 except for the case of n = 1. The latter behavior is in agreement with the results 
obtained in Ref. (l(| within the QCD sum rules. It is interesting to observe, that vanishing of <j>^ s (u) at the end 
points is correlated with the fact that (f>£ v (u) is concave at u ~ 1 and u ~ for n > 1. 

The paper is organized as follows: in Section [n] we shortly summarize the results of Refs. j5fl(| concerning the 
relation of the k T moments of the pion wave functions to the mixed quark-gluon condensates. In Sect. Ill, following 



Ref. J 13], we explain how to evaluate the loop integrals with the momentum dependent quark mass Mjk) given by 




In Sect. IV we construct the continuation prescription discussed above and calculate the gluon condensate. In 



we give the numerical results and Sect. VI contains summary and conclusions. Some of the results presented 
here have been partially reported in Ref. P6| . 
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II. PION WAVE FUNCTIONS AND THEIR RELATION TO THE VACUUM CONDENSATES 



We shall be dealing with the leading twist axial-vector (AV) and twist 3 pseudo-scalar (PS) wave functions defined 
as follows 

OO 

0f>) = _I_ / ^ e --( 2 -- 1 )(^)(0|^(r i TM75^(-nr)|7r + (P)>, 

l\ / 2F 7r J 7T 

— oo 

oo 

tf s (u) = _(„P) f — e--( 2 "- 1 )(" F )(0|^(nr)i75V(-nT)|7r+(P)). (4) 



V2 



qq) 



where we have chosen n = (1, 0, 0, —1) as a light-cone vector parallel to z M = rn p and n = (1, 0, 0, 1) parallel to P^. 
For any 4-vector v we have: 

+ - v+ ~ v ~ _l , r > 

v = n ■ v, v =n-v, = — + — + . (5) 

Both wave functions are normalized to 1. The normalization condition for <fi^ yields therefore the expression for 
(qq). Evaluation of the norm of (j>^ in the present model gives: 

, , T [ d * k n FTTTTTT, TvT k(k - P) - M(k)M(k - P) 

For P — > in the local case M{k) = M and in Euclidean metric Eq.(^) reduces to 

The last expression in (Q) is well defined only in Euclidean metric and Eq.(^) may serve as a prescription how it 
should be continued to the Minkowski space. In fact, as we shall explicitly see, (qq) of Eq.(^) does not depend on P. 
Our approach enables to calculate not only the distributions in u but in fact full wave functions ip v both in u and 

OO 

2 „/. („, ;„2 - 



<M«) = J dk T i(> w (u,k T ) (8) 



both for AV and PS channels. Integrating first over u gives the k\ distribution 

l 

07r(fer) = duipT,(u,k^). (9) 







It has been shown in Refs. [||^] that moments of w (fc^) are given in terms of the mixed quark-gluon condensates 

, 2 > 5 (zggq-Gg) , 2 . _ 1 (igqa-Gg) 

Here is a gluon field strength and 

c • G = a ^C" , G^ u = —G a ^ v . (11) 

Note that formula for (fey) ..y is an approximate one since the soft pion theorems used to derive ([!(]) apply strictly 
only in the pseudo scalar channel. Equations ( |To| ) predict 

m AV = ~ (4) PS . (12) 
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For (ftj-) we have |J: 



'Ay 



32 (V>V>) 72 (^) 



_ 1 (g^a-Ga-G^) l (g 2 ^) 

where for (fcr/^y the same reservations as for {k^) Av hold. Unfortunately equations (13) are almost linearly 
dependent. Indeed, coefficients in front of the condensates are of the order of 1/10, whereas the determinant of (13) 
is of the order 6 x 10~ 4 . Therefore it is difficult to disentangle the condensates from the knowledge of unless 
the latter are known with very high accuracy. 



III. PION WAVE FUNCTION WITH MOMENTUM DEPENDENT CONSTITUENT QUARK MASS 

In our previous paper [ p"3[ we have shown how to handle the loop integral, like the one in Eq.(^|), with momentum 
dependent mass M(k) given by Eq.(|). To this end one introduces the light cone integration variables 



d 4 k = -dk+dk-d 2 k T , 
2 



k 



+ _ 



iP+. 



(14) 



The key point of this analysis consists in finding the integration contour in the complex k plane. The poles of the 
fermion propagator are given by the zeros of the following function: 



M 



A 2 



which after pulling out the factor A 2 \k 2 /A 2 

'k 2 



1 + ie 

4rt+l 



k 2 - A 2 + ie 

in 



ie = 



reduces to 



A 2 



- 1 



A 2 



- 1 



in 



M 2 
A 2 " 



= 



and similarly for k — > k — P. Here 

k 2 = uP+k- 



and {k-Pf = -(l-u)P+r 



k'-p 



Eq.(|15[) should be understood as an equation for k 
Generally an equation of the form 



z 4n+l + z 4n _ ^2 







(15) 



(16) 



(17) 



with /j 2 = M 2 /A 2 has An + 1 complex solutions, which in the following will be denoted as where i = 1, . . . , An + 1. 
These solutions depend on the specific value of /x 2 and have to be calculated numerically. Then the poles of the quark 
propagators in k~ can be found from the following relations: 



fc 2 
A 2 

{P-kf 
A 2 



- 1 

- 1 



(18) 
(19) 



We have shown in Ref. |L3| that all the poles corresponding to Eq.(18) should lie below the integration contour whereas 
the poles corresponding to Eq.(19) above. As a result the dk~ integrals yield real wave functions, which vanish for u 
outside the region < u < 1. Moreover for A — > oo, i.e. for a constant M (k), this prescription reduces in a continuous 
way to the standard one of Feynman. 

With this prescription the calculations are rather straightforward and we obtain: 



^ V (u,k 2 ) 



1 N C M 2 ^ 



A 2 (2^) 2 F 2 .f- IlIk t + l + z t u + z k (l-u) : 



(20) 
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(where t = k 2 -,/ K 2 ). Factors /, 



N C M ^ z 3n z 3n (1 + a ±^ ) _^n 4 , 

(2tt) 2 (gg) /l/fc i + 1 + uz 4 + (1 - 



4n+l 



obey the following properties 



4n+l 



i=l 



for to < 4n 

1 for m = 4n 



which are crucial for the convergence of the d< integrals. The dt integration gives: 

N r M 2 



(2tt) 2 F 2 



£ [zfzfu + Z f"4(l - «)] In (1 + z lU + z k (l - u)) 



i,k 



and 



N C MA 2 



zf n zl n (l H - ) — ^ 2 z"z£ 



ln(l + uzi + (1 — 



(21) 



(22) 



(23) 



(24) 



(25) 



The normalization condition for <fr£ v is used to fix the value of the cutoff parameter A for given n with F„ fixed at the 
physical value of 93 MeV. Then the normalization condition for ip^ s can be used to calculate the quark condensate. 



IV. GLUON CONDENSATE 

The Euclidean formula for the gluon condensate in the instanton model of the QCD vacuum reads |l2| , p4| ]: 

a \ f d 4 k E M 2 (k E ) , x 

-GG) = 32N C / --^72 Sr, ' r . 26 

7T / J (2tt) 4 k% + M 2 {k E ) V ; 

In the instanton model, assuming that the gluon condensate is known, Eq. ( p6| ) should be understood as a gap equation 
for the value M = M(0), since the shape of the k dependence is uniquely determined as a Fourier transform of a 
fermion zero mode in the presence of the instanton configuration. Alternatively, for given M, one can use Eq.(|26|) to 
calculte the gluon condensate. 

A comparison of equations (0) and (|2^) for a constant constituent mass M (assuming some Euclidean cutoff) 
suggests the following relation between the quark and gluon condensates: 

M(qq) = -c(^GG). (27) 

Here M. = M (constituent quark mass) and c = 1/8. Let us remind that relation (|27]) follows from the instanton 
model of the QCD vacuum in the zero mode approximation and in the chiral limit (to = 0) [^4|. It would have to 



be substantially modified for heavy current quark mass to |27 2§[]. In Ref. |29j a similar relation has been obtained 
for an arbitrary self-dual gluonic field (in a dilute gas approximation) with, however, M = to (current quark mass of 
a light quark). In contrast to the previous case such a relation is not satisfied phenomenologically and moreover it 
cannot be true in the chiral limit. On the other hand an expansion for large current quark mass m q |]30|| yields in the 
leading order relation ( p^ ) with M. = m q and c = 1/12 pl|. It has been argued in Ref. |H| that this relation holds 
also for M. = M (constituent quark mass). Whether this identification can be theoretically justified is not obvious, 
since M being of the order of 300 MeV is probably not large enough to be considered as heavy. It is, however, beyond 
the scope of this paper to discuss here the differences between the approaches of Refs. ]2j|, |u) and Q (see, however 
[pU). We simply use Eq.(|26|) as derived in Ref. 0] and calculate the gluon condensate. 
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The first estimate of the gluon condensate given already in Ref. ||, (330 MeV) 4 , is nowadays believed to be too 
small. In fact different authors advocated different values for this quantity (see e.g. Fig. 4 of Ref. (3^ and references 
therein). For the purpose of this work we adopt more recent estimate from Refs. [ |34| 



-GG 



(393±|i MeV) 4 . 



(28) 



In our case we shall simply use Eq.jpq) to calculate and compare with Eq. (J28[) . 

It is tempting to use our experience from the wave function calculation to continue Eq. (^6|) to the Minkowski space. 
Indeed equations (^6|) and (0) differ only by one power of y/M(k)M(k — P) (apart from trivial numerical factors). 
This suggests the following generalization of (p6h: 



-GG) = i32N r 

7T 



M(k)M(k- P) 



k(k- P) - M(k)M(k- P) 



(29) 



(2tt) 4 *" v "'"' v ^ " ' (k 2 - M 2 (k)){(k - P) 2 - M 2 (k - P)) ■ 

For a constant M Eq.(|2^) can be easily obtained from ( ^6| ) by first continuing to the Minkowski metric and second 
by shifting the integration variable by a constant fourvector P (with P 2 = 0): 

la ^\ 00. , f d 4 k M 2 f d 4 k ( M 2 M 2 

-GG) = z327V c J (2?r)4 fc2 _ M2 = z!6iV c J ^ {j^W 2 + (k-Pf-M 2 



i32N c 



d 4 k M 2 (k(k — P) — M 2 ) 



(2tt) 4 (k 2 - M 2 )((fc - P) 2 - M 2 ) ' 

Accepting Eq.(^9|) as a continuation from the Euclidean to the Minkowski metric we can calculate the gluon condensate 
for the cutoff function of Eq. (pi) : 



-GG 



8N C M 2 A 2 
(2tt) 2 



du dt > fi fh - 

<^ JlJk t+l + u Zl + (l-u)z k 



Z. 



2n r ,2n 



(30) 



The dt integration is trivial (note that the large t part vanishes due to the properties of the sums of the fi factors) . 
Also the du integral is straightforward. Finally we obtain 



a \ _ 8N C M 2 A 2 
n) ~ (2tt)2 



i^k 



i 

Zfcs 2I f(l + Zi)\n{\ + Zi) - (1 + z fe )ln(l + z k ) 



- 1 



Zi - z k 



(31) 



P-k 




FIG. 1. Pion self-energy -iE^-(P) 



It is interesting to observe that the Minkowskian integral (|29j) is in fact proportional to the pion self-energy due to 
the quark loop calculated in the present effective model. Indeed 



2N r 

EJP) = -i—^r \ d 4 kM(k)M(k - P)Tr 
F n J 

which yields an interesting relation 



1 



1 



75 



:75 



jt. - M{k) 1 {jt-P)- M(k - P) 



-GG)=AF^(P). 



(32) 



(33) 



Let us note that in the nonlinear model which we are dealing with, there exists another loop diagram with two pions 
coupled to a fermion line at the same point, which cancels out ( |32"| ) so that the total pion self-energy is zero as it 
should be. 



G 



V. NUMERICAL RESULTS 



A. Quark and gluon condensates 



Our numerical results are presented in Table || For 4 different values of the constituent mass M = M(0) and 
different n (see Eq.([|)) we have adjusted the cutoff parameter A by imposing the normalization condition on <p^ v ■ 
In fact, as discussed at length in Ref. |l3| the leading twist pion wave function <f>^ v (u) does not change any more if 
we increase n above 5. On the other hand for n > 5 the cutoff function (^), if continued to the Euclidean metric, 
starts to deviate significantly from the one obtained in the instanton model. Therefore we have chosen to work with 
n max = 5. From our previous study it seems that the best agreement with the recent analysis of the CLEO data is 
obtained for M = 325 MeV and n = 2 - 5 or M = 350 MeV and n = 2. 

In the last two columns of Table | we have displayed the values of the quark and gluon condensates. The gluon 
condensate is quite insensitive to the parameters of the cutoff function and also to the value of the constituent quark 
mass M . The numerical value being of the order (380) 4 -(400) 4 MeV 4 agrees perfectly with the phenomenological 
value. In contrast, the quark condensate is sensitive to the value of n and the numerical value lies between — (258) 3 and 
— (338) 3 MeV 3 overshoots the phenomenological value of — (250) 3 MeV 3 . The reason for this sensitivity is relatively 



broad k\ distribution of (f)^ s (k^). We shall come back to this point in Sect.VC 



B. Pseudo-scalar pion wave function 



The function <p^ s (u) has been calculated within the QCD sum rules in Refs. |L5 16 1. It had a u— shape and did not 
vanish at the end points. In Fig. 2 we plot our results for <p^ s (u). It can be seen, that for n = 1 <p^ s (1) — <fi^ s (0) ^ 0, 
whereas for n > 2 <fi^ (1) = cf)^ (0) = 0. The end point behavior is governed by the sum 



6^(1) ~5>ln(l + 



Jin Jin , 1 3n+l in . 1 in in+1 ,,2 ran 
z i z k ' 2 i fe + 2 k 



(34) 



which vanishes due to the property (123) except for n = 1 where 3n + 1 = An. 



TABLE I. Quark and gluon condensates for different choices of cutoff parameters M(0), n and A. 



M(0) 
Mev 


n 


A 
MeV 


MeV 3 


<^GG) 
MeV 4 


325 


1 

2 
5 


1249 
1862 
3033 


(328) d 
(294) 3 
(280) 3 


(403)" 
(396) 4 
(393) 4 


350 


1 
2 
5 


1156 
1727 
2819 


(318) 3 
(284) 3 
(271) 3 


(399) 4 
(392) 4 
(389) 4 


375 


1 
2 
5 


1081 
1621 
2649 


(309) 3 
(277) 3 
(264) 3 


(395) 4 
(388) 4 
(386) 4 


400 


1 
2 
5 


1020 
1543 
2512 


(303) 3 
(271) 3 
(258) 3 


(392) 4 
(386) 4 
(383) 4 
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M = 325 MeV 



M = 350 MeV 





FIG. 2. Pseudoscalar pion wave function <^ s (w) for M = 325, 350, 375 and 400 MeV, for n = 1 (solid), 2 (long dashed) and 
5 (dashed). 



C. fcy dependence and mixed vacuum condensates 

The large fc|> asymptotics of the wave functions and of the appropriately defined k\ distribution of the gluon 
condensate 

a GG) = / dk 2 T G(k 2 T ) (35) 



7T 







is as follows 



^ \ 4n+l / -, \ 2n / 1 \ 4n 



#r(#)~(^-J > 7T - <5(4)~ t^) (36) 



Therefore for example for n = 1 already the first moment of 0^ >s is divergent. 

As a kind of remedy for this divergence let us try to introduce a cutoff for dk 2 -, integration T — K^/A 2 . In Table 
II we list k\ moments, and (— GG) for T chosen in such a way that the quark condensate equals (qq) — —(250 
MeV) 3 . Interestingly, by fixing (qq) with one parameter T we are able to reproduce ratio (|l^) with high accuracy 
for almost all sets of parameters. Unfortunately the ratio {k?) AV I \^t/ps IS ^ ess s t am e: 0.27 — 0.43. Moreover it is 
more than two times smaller than the one obtained by means of Eqs.(13) in Ref. Q. 

One should, however, keep in mind that th e absolute values of k\ moments, especially in the PS channel, depend 



strongly on T. This is illustrated in Table |I| for M = 350 MeV and n = 2. 

One more remark is here in order. We have not readjusted A since F n (T — 0.25) ~ 92 MeV instead of 93, a negligible 

change. In Fig.4 we plot T dependence of the ratios F W (T)/F n , {/(qq) (T)/^qj and {J (%GG) (T)/^/(fGG) 
for n = 2. We see that while F n and gluon condensate saturate relatively fast, the quark condensate (that is 
J dk^<j)^ s (k 2 -,)) saturates very slowly. 
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M = 350GtV 



■e- 




kl [GeV ] 



FIG. 3. Dependence on h\ of fa b {u) for M = 350 MeV, for n = 1 (solid), 2 (long dashed) and 5 (dashed) and of <££ v (i 
for all n = 1 ... 5 (dot-dashed) 



TABLE II. gluon condensate and fcy moments of axial- vector and pseudo-scalar pion wave functions with transverse 
integration cutoff fixed by requirement that (qq) = —(250 MeV) 3 . 



M 


n 


T 


K T 






y/(k*) 


(4) 4V 

( fc x) PS 


V<4> 




MeV 






MeV 


MeV 


MeV 


MeV 

AV 


PS 


MeV 

AV 


PS 




1 


0.374 


764 


89.5 


378 


356 


476 


0.56 


428 


529 


0.43 


325 


2 


0.191 


813 


90.7 


381 


370 


495 


0.56 


445 


554 


0.42 




5 


0.079 


852 


91.4 


382 


378 


508 


0.56 


456 


570 


0.41 




1 


0.47 


792 


90 


380 


363 


487 


0.55 


436 


544 


0.41 


350 


2 


0.25 


864 


91.6 


382 


378 


513 


0.54 


455 


577 


0.39 




5 


0.105 


913 


92 


383 


385 


527 


0.53 


465 


597 


0.37 




1 


0.58 


823 


91 


380 


368 


499 


0.54 


442 


558 


0.39 


375 


2 


0.31 


902 


92 


382 


381 


524 


0.53 


459 


592 


0.36 




5 


0.14 


991 


92.6 


383 


388 


548 


0.50 


471 


626 


0.32 




1 


0.69 


847 


91.5 


380 


371 


507 


0.53 


445 


569 


0.37 


400 


2 


0.39 


964 


92.5 


382 


383 


540 


0.50 


464 


615 


0.32 




5 


0.189 


1092 


92.9 


382 


389 


569 


0.47 


474 


659 


0.27 



TABLE III. Same as in Table [ll| for M = 350 MeV and n = 2 compared with the results with no fcy cutoff. 



t f w tfiioG) VW) VW) 

MeV MeV MeV MeV ' MeV 

AV PS AV^ PS_ 

0.25 250 91.6 382 378 513 0.54 455 577 

oo 284 93 392 420 921 0.21 540 1360 







1,0-1 



0,9- I' 

I / 
J / 

0,8- | / 
i 

0,7- \ j 

0,6 -i / 

I j 
' : i 

0,5 -F 1 1 1 1 1 1 . 1 1 1 

0,0 0,2 0,4 0,6 0,8 1,0 



FIG. 4. Dependence on the cutoff T of the following ratios: F vr (T)/F 7r (solid), {/ ' (qq) (T)/y/(qq) (dashed) and 
(f GG) (T) / d ( f GG) (dash-dotted) for n = 2. 

In order to estimate the value of the mixed condensate of dimension 5, (igqcr ■ Gq), we choose to use the first 
equation of Eqs.(|To|). Although this equation is only an approximate one, in the case of no transverse momentum 
cutoff it gives relatively stable results due to the fast convergence of the k\ integrals. In the case of the finite transverse 
momentum cutoff any of the two equations ([if]) can be used, since the relation ([[2]) is fulfilled with good accuracy. 
In the case of no cutoff (i.e. T ~ oo) we get that (for (qq) from Table [j]): 

(igqcr ■ G q) = - 450 5 to - 550 5 MeV 5 . (37) 

The value of ( |37| ) decreases with increasing M and with increasing n. Interestingly, for the parameters which are 
closest to the original instanton model, M — 350 MeV and n — 2, we get —(493 MeV) 5 in perfect agreement with 
the direct calculation of (igqcr ■ Gq) in the instanton model |n|, which gives —(490 MeV) 5 . In the case of the finite 
cutoff the value of (igqcr ■ Gq) is smaller and almost independent of the specific set of parameters (^): 

(igqa ■ Gq) = - 427 5 to - 443 5 MeV 5 . (38) 

This range of values is compatible with the result obtained within the QCD sum rules [^| at low normalization point: 
-(416 MeV) 5 . 



VI. SUMMARY AND CONCLUSIONS 

In this paper, following our previous work and Refs. |l7j , |l8| ] , we have studied a simple, nonperturbative model 
in which quarks acquire a momentum dependent dynamical mass (|J). In Ref. |l3| we have employed this model to 
calculate the axial-vector leading twist pion light cone wave function, whereas here we have extended our calculations 
to the case of the pseudo-scalar, twist 3 wave function (f>^ s , which is normalized to the quark condensate. Since the 
model parameters are fixed by the normalization of the axial-vector wave function we could use the normalization 
condition for (f)^ s to calculate the quark condensate. One of the advantages of the present model is that one can 
calculate not only the longitudinal quark distribution in the pion, but also k\ distributions both in the AV and PS 
channels. Therefore we could use the relations between k\ moments and mixed quark gluon condensates j^J(| to 
estimate the latter in the present model. 

One remark is here in order. The advantage of our approach is its simplicity, both conceptual and technical, 
which makes it possible to compute many quantities analytically. One might therefore suspect that the model is 
oversimplified and that our assumptions are too crude to be in general true. For example the model has no trace of 
confinement. It is therefore of importance to perform various tests in order to gain confidence in the model as well 
as to find its limitations. This study shows that in fact the model works much better than one might have initially 
expected. As shown in Ref. |l3) 4>^ v is compatible with recent CLEO data ]3q] ; quark condensate calculated here 
from the normalization condition of (f>^ s comes out reasonable, although bigger than the phenomenological value. The 
gluon condensate, which we have computed from the gap equation (Um guessing the continuation to the Minkowski 
space, comes out very well. The mixed quark-gluon condensate of dimension 5 calculated from the value of (fcf-) Av 
(see Eq. (fL0|) ) is in a surprising agreement with the direct evaluation of (igqcr ■ G q) in the instanton model |To[ . 
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An obvious limitation of the model is due to the k T asymptotics of 4>^ s and also <f>£ v (|36|). This powerlike behavior 
makes higher k T moments divergent, and therefore higher mixed condensates uncalculablc. For (fr£ v the power is still 
relatively high, whereas for (j>^ s it is rather low. This is the reason why the quark condensate comes out higher that 
the phcnomcnological value and depends (within a factor of 2) on model parameters. Since there exists another wave 
function normalized to (qq), namely 0^ Jig], it would be interesting to calculate (qq) from (f>° and compare with the 
present results p6[ . As a remedy for the broadness of (j>^ s (k T ) we have employed the k T cutoff, which does not affect 
much the gluon condensate and (i.e. 4>^ v ) but allows to reproduce the phenomenological value of (qq) = —(250 
MeV) 3 . Then the ratio (k T ) Av I (&t) ps ~ 5/9 is in agreement with the analysis of Ref. j5j. Unfortunately one single 
cutoff is probably not enough to make all higher moments reliable. 

Both quark and gluon condensate depend on the scale Qq. In model calculations it is, however, hard to define 
precisely what value should be taken for Qq. In the instanton model of the QCD vacuum the scale is determined 
by the average inverse size of the instanton l/jo = 600 MeV. It is therefore natural to expect that also in our case 
the normalization scale is of the same order. Its connection to the scale A entering Eq.(j^) is, however, by no means 
straightforward since in fact we use various values of A and n to approximate the same function M(k) evaluated in 
the instanton model for 1/p = 600 MeV. One way to determine Qq in the present approach is to identify Qo = Kt, 
where Kt is the transverse cutoff displayed in Table [nj. Indeed, this is precisely the way how the normalization scale 
is defined for the light cone wave functions (see e.g. Ref. Numerically it comes out that 760 < Qo < 1100 MeV. 

The twist 3 pion wave function <p^ s , which we have plotted in Fig. 2 is theoretically an interesting object, since it 
tends asymptotically to 1, and therefore might be enhanced in certain kinematical regions. The calculations of <fi^ s 
from the QCD sum rules suggest that it does not vanish at the end points |l5|;[l6|]. We find this kind of behavior only 
for n = 1 (see Eq.(^)), whereas for n > 1 (f>^ s vanishes at the end points. Interestingly, the vanishing of (f>^ s for 
u = 0, 1 is correlated with the nonconvexity of the axial-vector wave function at the end points, which for u — > (or 
1) behaves like u n (or (1 — u) n ). It would be worthwhile to check if this correlation holds also in other models. The 
fact, that the axial- vector pion distribution amplitude may be concave at the end points has been already pointed out 



in Ref. |36|] and confirmed by model calculations within the framework of the nonlocal sum rules 37 
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